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A brief review of the stability of rotating relativistic stars is followed by a more detailed 
discussion of recent work on an instability of r-modes, modes of rotating stars that have axial 
parity in the slow-rotation limit. These modes may dominate the spin-down of neutron stars 
that are rapidly rotating at birth, and the gravitational waves they emit may be detectable. 

§1. Introduction 

A series of recent surprises appear dramatically to have improved the likelihood 
that the spin of rapidly rotating, newly formed neutron stars is limited by a nonax- 
isymmetric instability driven by gravitational waves - and that the emitted waves 
may be detectable. 

The first of these was the discovery that the r-modes, rotationally restored modes 
that have axial parity for spherical models, are unstable in perfect fluid models with 
arbitrarily slow rotation. First indicated in numerical work by Andersson^^, the 
instability is implied in a nearly newtonian context by the newtonian expression 
for the r-mode frequency (3-1), and a computation by Friedman and Morsink^^^ 
of the canonical energy of initial data showed (independent of assumptions on the 
existence of discrete modes) that the instability is a generic feature of axial-parity 
perturbations of relativistic stars. 

Studies of the viscous and radiative timescales associated with the r-modes 
(Lindblom et al. Owen et al. ^'^\ Andersson et al. Kokkotas and Stergioulas ■^^-^ , 
Lindblom et al. '^'')) have revealed a second surprising result: The growth time of 
r-modes driven by current-multipole gravitational radiation is significantly shorter 
than had been expected, so short, in fact, that the instability to gravitational radi- 
ation reaction easily dominates viscous damping in hot, newly formed neutron stars 
(see Fig. 1 below). As a result, a neutron star that is rapidly rotating at birth 
now appears likely to spin down by radiating most of its angular momentum in 
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gravitational waves. (See, however, the caveats below.) 

Nearly simultaneous with these theoretical surprises was the discovery by Mar- 
shall ct.al. of a fast (16ms) pulsar in a supernova remnant (N157B) in the Large 
Magellanic Cloud. Estimates of the initial period put it in the 6-9ms range, implying 
the existence of a class of neutron stars that are rapidly rotating at birth. Hence, 
the newly discovered instability appears to set the upper limit on the spin of the 
newly discovered class of neutron stars. 

We discuss these developments and the uncertainties that accompany them be- 
low, in the context of a review of the gravitational- wave driven instability of rotating 
stars. The structure and stability of rotating relativistic stars has recently been re- 
viewed in detail (Stergioulas Friedman Friedman and Ipser-*^^)), and a 
general discussion of the small oscillations of relativistic stars may be found in a 
recent review article by Kokkotas^^^ (see also Kokkotas and Schmidt ■^^-'). 

The excitement over the r-mode instability has generated a large literature in the 
past two years. (Andersson -"^^ Friedman and Morsink^^^ Kojima^^\ Lindblom et 
al.^^\ Owen et al.^^\ Andersson, Kokkotas and Schutz^), Kokkotas and Stergioulas 

Andersson, Kokkotas and Stergioulas Madsen^'^), Hiscock^^), Lindblom and 
Ipser , Bildsten , Levin ^'^^ , Ferrari et al. , Spruit , Brady and Creighton , 
Lockitch and Friedman Lindblom et al. ^"^^ Beyer and Kokkotas Kojima and 
Hosonuma^^), Lindblom -^^^ , Schneider et al. RezzoUa et al. Yoshida and Lee 
) [Earlier studies of the axial-parity oscillations of models of the neutron star crust 
were reported by van Horn^-'^) and by Schumaker and Thorne^^^); and Chandrasekhar 
and Ferrari discussed the resonant scattering of axial wave modes and the coupling 
between axial and polar modes induced by stellar rotation .] Although we mention 
many of the recent contributions listed here, our summary is too short to allow much 
detail. 

§2. Gravitational- Wave Driven Instability of Rotating Relativistic Stars 

Only recently have the oscillation modes of rotating relativistic stars begun to 
be accessible to numerical study (see below). Early work on the perturbations of 
such stars focused mainly on the criteria for their stability, and led to the discovery 
that all rotating perfect fluid stars are subject to a nonaxisymmetric instability 
driven by gravitational radiation. The instability was found by Chandrasekhar 
for the I = m = 2 polar mode of the uniform-density, uniformly rotating Maclaurin 
spheroids. Although this mode is unstable only for rapidly rotating models, by 
looking at the canonical energy of initial data with arbitrary values of m, Friedman 
and Schutz^'') and Friedman showed that the instability is a generic feature of 
rotating perfect fluid stars, that even slowly rotating perfect-fluid models are formally 
unstable. 

For a normal mode of the form e*(^*+"*'^) this nonaxisymmetric instability acts 
in the following manner: In a non-rotating star, gravitational radiation removes 
positive angular momentum from a forward moving mode and negative angular 
momentum from a backward moving mode, thereby damping all time-dependent, 
non-axisymmetric modes. In a star rotating sufficiently fast, however, a backward 



Gravitational- Wave Driven Instability 



3 



moving mode can be dragged forward as seen by an inertial observer; and it will then 
radiate positive angular momentum. The angular momentum of the mode, however, 
remains negative, because the perturbed star has lower total angular momentum 
than the unperturbed star. As positive angular momentum is removed from a mode 
with negative angular momentum, the angular momentum of the mode becomes in- 
creasingly negative, implying that its amplitude increases: The mode is driven by 
gravitational radiation. 

The conclusion, that a mode is unstable if it is prograde relative to infinity and 
retrograde relative to the star is equivalent to requiring that its frequency satisfies 
the condition, 

a{a + mf2)<0. (2-1) 

For the polar f- and p-modes, the frequency is large and approximately real. Condi- 
tion (2-1) will be met only if \mf2\ is of order \a\, so that for a given angular velocity 
the instability will set in first through modes with large m. 

The instability spins a star down by allowing it to radiate away its angular 
momentum in gravitational waves. However, to determine whether this mechanism 
may be responsible for limiting the rotation rates of actual neutron stars, one must 
also consider the effects of viscous damping on the perturbations. Dctweiler and 
Lindblom -"^^^ suggested that viscosity would stabilize any mode whose growth time 
was longer than the viscous damping time, and this was confirmed by Lindblom and 
Hiscock . Recent work has indicated that the gravitational- wave-driven instability 
can only limit the rotation rate of hot neutron stars, with temperatures above the 
superfluid transition point, T ~ lO^K, but below the temperature at which bulk 
viscosity apparently damps all modes, T ~ lO^^K (see Fig. 1). (Ipser and Lind- 
blom Lindblom and Lindblom and Mendell^^^) Because of uncertainties in 
the temperature of the superfluid phase transition and in our understanding of the 
dominant mechanisms for effective viscosity, even this brief temperature window is 
not guaranteed. 

The instability appears to arise in modes with both polar and axial parity, 
with the distinction defined as follows. The spherical symmetry of a nonrotating 
star implies that its perturbations can be divided into two classes, polar and axial, 
according to their behavior under parity. Where polar tensor fields on a 2-sphere 
can be constructed from the scalars YJ"^ and their gradients VY^"^ (and the metric 
on a 2-sphere), axial fields involve the pseudo- vector f x Vy^™", and their behavior 
under parity is opposite to that of Yj^. That is, axial perturbations of odd I are 
invariant under parity, and axial perturbations with even I change sign. Because a 
rotating star is also invariant under parity, its perturbations can also be classified 
according to their behavior under parity. Although I is well-defined only for modes 
of the spherical configuration, if a mode varies continuously along a sequence of 
equilibrium configurations that starts with a spherical star and continues along a 
path of increasing rotation, the mode can be called axial (polar) if it is axial (polar) 
for the spherical star. 

It is useful to subdivide stellar perturbations according to the physics dominat- 
ing their behaviour. This classification was first developed by Cowling for the 
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Temperature (K) 

Fig. 1. Critical angular velocity vs. temperature for an n = 1.0 polytrope. Above the solid curve, 
the star rotates rapidly enough for its fastest growing (Z = m = 2) r-mode to be unstable, 
whereas below the curve all modes are damped by viscosity. The dashed curve does not include 
the effects of bulk viscosity, while the solid curve does; however, for more accurate calulations of 
the bulk viscosity contribution see Lindblom, Mendell and Owen and Andersson, Kokkotas 
and Schutz . The dot-dashed curve shows the evolution of a rapidly rotating neutron star as it 
cools and spins down due to the emission of gravitational waves. [This figure is due to Lindblom, 
Owen and Morsink and is reproduced here by permission of the authors.] 



polar perturbations of newtonian polytropic models. The p-modes are polar-parity 
modes having pressure as their dominant restoring force *^ They typically have large 

pressure and density perturbations and high frequencies (higher than a few kilohertz 
for neutron stars). The other class of polar-parity modes are the g- modes, which 
are chiefly restored by gravity. They typically have very small pressure and density 
perturbations and low frequencies. Indeed, for isentropic stars, which are marginally 
stable to convection, the g-modes arc all zero- frequency and have vanishing perturbed 
pressure and density. Similarly, all axial-parity perturbations of newtonian perfect 
fluid models have zero frequency in a non-rotating star. The perturbed pressure and 
density as well as the radial component of the fluid velocity are all rotational scalars 
and must have polar parity. Thus, the axial perturbations of a spherical star are 
simply stationary horizontal fluid currents. 

The analogues of these modes in relativistic models of neutron stars have been 
studied by many authors. More recently, an additional class of outgoing modes has 
been identified that exist only in relativistic stars. Like the modes of black holes. 



The lowest p-mode for each value of I and m is termed an f-mode or fundamental mode. 
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these are essentially associated with the dynamical spacetime geometry and have 
been termed w-modes, or gravitational wave modes. Their existence was first argued 
by Kokkotas and Schutz^^\ The polar w-modes were first found by Kojima^^^ as 
rapidly damped modes of weakly rclativistic models, while the axial w-modes were 
first studied by Chandrasekhar and Ferrari^'') as scattering resonances of highly 
relativistic models. (See the reviews by Kokkotas and Kokkotas and Schmidt ^^).) 
All w-modes appear to be stable, and if they are analogous to the modes of non- 
interacting zero-rest-mass test fields on the geometry of a rotating star, they should 
be unstable only when an ergosphere is present; but we are not aware of a careful 
study to rule out w-mode instability. 

Until recently, the polar p-modes were expected to dominate the CFS instability 
through their coupling to mass multipole radiation. The relativistic stability points 
of p-modes have recently been found by Stergioulas and Friedman ^^^'^^^ and, in 
the Cowling approximation, by Yoshida and Eriguchi^^). This last reference also 
provides the first computation of the modes of rapidly rotating relativistic stars. 
The relativistic computation shows an instability of each at significantly smaller 
values of dimensionless measures of rotation than is seen in nearly newtonian stars. 
Particularly striking is the fact that the l=m=2 bar mode is unstable for relativistic 
polytropcs of index n=1.0. The classical Newtonian result for the onset of the bar 
mode instability, n < 0.808) is replaced by ra < 1.3 for relativistic polytropes. 

Fig. 2 shows the critical angular velocity above which the lowest / = m p-modes 
are unstable in perfect fluid models. For most realistic equations of state, Morsink, 
Stergioulas and Blattnig^^) find that the I = m = 2 bar mode is unstable in IAMq 
neutron stars for angular velocities above O.SOk to 0.95^?;^. 

The competition between viscosity and gravitational radiation is described by 
the dissipation equation, 

111 1 , , 

- = + + , (2-2) 

T T^GK '''shearviscosity "njulkviscosity 

with r the e-folding time for each process. (The analysis sketched below can be found 
in Lindblom ct al'^^^; see also Ipser and Lindblom^^)). When the energy radiated 
per cycle is small compared to the energy of the mode, the imaginary part of the 
mode frequency is accurately approximated by the expression 

r" 2Edt' ^^-^^ 
where E is the energy of the mode as measured in the rotating frame, 



(fx. (2-4) 



We have. 



dt 



l>2 



{2r,6a'''dal, + CdeSe*), (2-5) 
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N=1.0 




Fig. 2. Critical ratio of angular velocity to Kcplcr velocity vs. a dimonsionloss central energy 
density ec for the m = 2, 3, 4 and 5 neutral modes of n = 1.0 polytropes. The largest value of tc 
shown corresponds to the most relativistic stable configurations, while the lowest €"c corresponds 
to less relativistic configurations. The filled circles on the vertical axis represent the Newtonian 
limit. 

where the dissipation due to gravitational radiation has coupUng constant 

47rG + + 2) 
' c2'+i;(/-l)[(2/ + l)!!]2' 

5(Jab and 59 are the coefficients of shear and bulk viscosity; and the corresponding 
coefficients r] and C, are estimated (Cutler and Lindblom^^)^ Sawyer '^^)) by 

and 
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Polar and axial radiation arise, respectively, from mass and current multipoles, 
Dim and Jim, given by the equations, 

Dim = J dVr^pYCm Jim = J dVr^pw ■ r x VY{^ (2-9) 

The additional factor of v in the current multipoles implies an additional factor 
of in the radiated energy of axial modes, and hence a smaller expected rate of 
radiation for the same multipole: For a mode of amplitude A = (displacement of 
fluid element)/-R, with R the stellar radius, we have 

^ ~ A^M^i?^; 2i+2 dE ^ ^2^2^2;+2 a+4_ (2-10) 

dt dt ^ ' 

The extra factor of and the fact that a is proportional to for slowly rotating 
stars led to the incorrect expectation that r-modes could be neglected. 

§3. The r-mode Instability and its Implications 

The unexpected apparent dominance of the r-modes arises from the fact that, 

where the polar I = 2 mode is unstable, if at all, only for Q ^ Q^, the I = m = 
2 r-mode is unstable (for perfect-fluid models) for arbitrarily small Q, The 
instability in models of slowly rotating, nearly newtonian stars, follows from the fact 
that the frequencies, 

satisfy the criterion (2-1), 

<^ + mQ) = < 0. (3-2) 

[A computation by Friedman and Morsink of the canonical energy of initial data 
showed (independent of assumptions on the existence of discrete modes) that the 
instability is a generic feature of axial-parity fluid perturbations of relativistic stars.] 

Not only are p-modes with low values of m unstable only for large values of 
f2, their frequencies are near zero, because the frequency of a mode vanishes at its 
instability point, and a star cannot rotate fast enough for the lowest p-modes to be 
far from their instability points. Thus where p-modes may limit the rotation of a 
newly formed star to 0.9f2K,^^^ *^ the r-modes are likely to limit rotation to less 
than 0.2ilK- Here is the maximum angular velocity of a uniformly rotating star, 
equal to the angular velocity of a satellite in orbit at the star's equator. 

The current picture that has emerged of the spin-down of a hot, newly formed 
neutron star can be readily understood in terms of a model of the r-mode instability 
due to Owen, Lindblom, Cutler, Schutz, Vecchio and Andersson (hereafter OLCSVA) 
Since one particular mode (with spherical harmonic indices I = m = 2 and 

*^ this is likely to decrease, perhaps to 0.8 Qk when the correct general-relativistic computation 
with viscosity is completed. 
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frequency a = — 4i7/3) is found to dominate the r-mode instability, the perturbed 
star is treated as a simple system with two degrees of freedom: the uniform angular 
velocity f2 of the equilibrium star, and the (dimcnsionlcss) amplitude a of the I = 
m = 2 r-mode. Initially, the neutron star forms with a temperature large enough 
for bulk viscosity to damp any unstable modes, T ^ lO^'^K; the star is assumed to 
be rotating close to its maximum (Kepler) velocity, fix ~ \/M/R-^. The star then 
cools by neutrino emission at a rate given by a standard power law cooling formula 
(Shapiro and Teukolsky ) . Once it reaches the temperature window at which the 
I = m = 2 r-mode can go unstable, the system is assumed to evolve in three stages. 

First, the amplitude of the r-mode undergoes rapid exponential growth from 
some arbitrary tiny magnitude. Using conservation of energy, (2-5) and the related 
equation for the rate of angular momentum loss to gravitational waves, 

^ = '^22' (3-3) 

with both E and J proportional to a^, OLCSVA derive the following equations for 
the evolution of the system in this stage. After expressing E and J in terms of Q 
and a 

df2 _ 2Q a^Q 

\^'^) 



dt ry 1 + a'^Q 
da a a 1 — c?Q 



(3-5) 



\tgr\ Tyl + a'^Q 

Here, tgr and Ty are, respectively, the timescales for the growth of the mode by 
gravitational radiation reaction and the damping of the mode by viscosity (see Sect. 
2.6). (The parameter Q is a constant of order 0.1 related to the initial angular mo- 
mentum and moment of inertia of the equilibrium star.) Since the initial amplitude 
a of the mode is so small, the angular momentum changes very little at first (Eq. 
(3-4)). That this stage is characterized by the rapid exponential growth of a is the 
statement that the first term in Eq. (3-5) (the radiation reaction term) dominates 
over the second (viscous damping). 

Eventually the mode will grow to a size at which linear perturbation theory is 
insufficient to describe its behavior. It is expected that a non-linear saturation will 
occur, halting the growth of the mode at some amplitude of order unity, although the 
details of these non-linear effects are poorly understood at present. When this satu- 
ration occurs, the system enters a second evolutionary stage during which the mode 
amplitude remains essentially unchanged and the angular momentum of the star is 
radiated away. During this stage OLCSVA evolve their model system according to 
the equations 

= K (3-6) 

df2 _ 2Q kQ 

^ ~ ~\tgr\1-i^Q ' 
where n is constant of order unity parameterizing the uncertainty in the degree of 
non-linear saturation. The star spins down by Eq. (3-7), radiating away most of its 
angular momentum while continuing to cool gradually. 
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When its temperature and angular velocity are low enough that viscosity again 
dominates the gravitational-wave-driven instability, the mode will be damped. Dur- 
ing this third stage, OLCSVA return to Eqs. (3-4)- (3- 5) to continue the evolution of 
their system. That the mode amplitude decays is the statement that the second term 
in Eq. (3-5) (the viscous damping term) dominates the first (radiation reaction), at 
this temperature and angular velocity. 

This three-stage evolutionary process leaves the newly formed neutron star with 
an angular velocity small compared with Qk- This final angular velocity appears 
to be fairly insensitive to the initial amplitude of the mode and to its degree of 
non-linear saturation. A final period P ^ 5 — 10ms apparently rules out accretion- 
induced collapse of white dwarfs as a mechanism for the formation of millisecond 
pulsars with P ^ 3ms. 

The r-mode instability has also revived interest in the Wagoner mechanism, 
in which old neutron stars are spun up by accretion until the angular momentum loss 
in gravitational radiation balances the accretion torque. Bildsten^^ and Andersson, 
Kokkotas and Stergioulas have proposed that the r-mode instability might succeed 
in this regard where the instability to polar modes seems to fail. However, the 
mechanism appears to be highly sensitive to the temperature dependence of viscous 
damping. Levin ^'^^ has argued that if the r-mode damping is a decreasing function of 
temperature (at the temperatures expected for accreting neutron stars, T ~ lO^K) 
then viscous reheating of the unstable neutron star could drive the system away from 
the Wagoner equilibrium state. Instead, the star would follow the cyclic evolution 
pattern depicted in Fig. 3. Initially, the runaway reheating would drive the star 
further into the r-mode instability regime (B-C in Fig. 3) and spin it down to a 
fraction of its angular velocity (C-D). Once it has slowed to the point at which the 
r-modes become damped, it would again slowly cool (D-A) and begin to spin up by 
accretion (A-B). Eventually, it would again reach the critical angular velocity for the 
onset of instability and repeat the cycle. Since the radiation spin-down time is of 
order 1 year, while the accretion spin-up time is of order 10^ years, the star spends 
only a small fraction of the cycle emitting gravitational waves via the unstable r- 
modes. This would significantly reduce the likelihood that detectable gravitational 
radiation is produced by such sources. On the other hand, if the r-mode damping is 
independent of - or increases with - temperature (at T ~ lO^K) then the Wagoner 
equilibrium state may be allowed (Levin ) . Work is currently in progress (Lindblom 
and Mendell ) to investigate the r-mode damping by mutual friction in superfluid 
neutron stars, which was the dominant viscous mechanism responsible for ruling out 
the Wagoner scenario in the first place (Lindblom and Mendell ) . 

Other uncertainties in the scenarios described above are still to be investigated. 
There is substantial uncertainty in the cooling rate of neutron stars, with rapid 
cooling expected if stars have a quaxk interior or core, or a kaon or pion condensate. 
Madsen^°) suggests that an observation of a young neutron star with a rotation 
period below 5 — 10ms would be evidence for a quark interior; but even without 
rapid cooling, the uncertainty in the superfiuid transition temperature may allow a 
superfluid to form at about 10^°K, possibly killing the instability. 

We noted above the expectation that the growth of the unstable r-modes will 



10 



J. L. Friedman and K. H. Lockitch 




Fig. 3. Levin's ^ cyclic evolution scenario for an accreting neutron star when r-rnodc damping is a 
decreasing function of temperature. The solid curve K-L-M shows the critical angular velocity 
vs. temperature, while the dashed curve A-B-C-D-A shows the evolution of the star. [This 
figure is due to Levin and is reproduced here by permission of the author.] 



saturate at an amplitude of order unity due to non-linear effects (such as mode-mode 
couplings); however, this limiting amplitude is not yet known with any certainty and 
could be much smaller. In particular, it has been suggested that the non-linear 
evolution of the r-modes will wind up the magnetic field of a neutron star, draining 
energy away from the mode and eventually suppressing the unstable modes entirely 
(RezzoUa, Lamb and Shapiro see also Spruit ^^•*). 

If the r-modes saturate near unity, a study by Brady and Creighton^^ follow- 
ing the work by Owen et al^^\ finds that newly formed neutron stars should be 
detectable by LIGO II with narrow banding out to about 8 Mpc, with uncertainty 
allowing a range of perhaps 4-20 Mpc. The Virgo cluster is then likely to be out 
of reach, and even with this optimistic assumption about the saturation amplitude, 
r-modes are most likely not to be detected until sensitivity passes the limits of ad- 
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vanced LIGO. 

Despite the recent intense interest in r-modes, they are not yet well-understood 
for stellar models appropriate to neutron stars. A neutron star is accurately described 
by a perfect fluid model in which both the equilibrium and perturbed configurations 
obey the same one-parameter equation of state. We call such models isentropic, 
because isentropic models and their adiabatic perturbations obey the same one- 
paramctcr equation of state. 

For stars with more general equations of state, the r-modcs describe the dynami- 
cal evolution of initial perturbations that have axial parity. This is not, however, the 
case for isentropic models. Early work on the r-modes focused on newtonian mod- 
els with general equations of state (Papalouizou and Pringle^^\ Provost et al.^^^, 
Saio^^\ Smeyers and Martens ) and mentioned only in passing the isentropic case. 
In isentropic newtonian stars, one finds that the only purely axial modes allowed are 
the r-modes with I = m and simplest radial behavior. (Provost et al.^^^*)) It is 
these r-modes only that have been studied (and found to be physically interesting) 
in connection with the gravitational-wave driven instability. 

The disappearance of the purely axial modes with I > m occurs for the follow- 
ing reason. (Lockitch and Friedman ^^^) Axial perturbations of a spherical star are 
time-independent convective currents with vanishing perturbed pressure and den- 
sity. In spherical isentropic stars, stars for which both star and perturbation are 
governed by a single one-parameter equation of state, the gravitational restoring 
forces that give rise to the g-modes vanish and they, too, become time-independent 
convective currents with vanishing perturbed pressure and density. Thus, the space 
of zero frequency modes, which generally consists only of the axial r-modes, expands 
for spherical isentropic stars to include the polar g-modes. This large degenerate 
subspace of zero-frequency modes is split by rotation to zeroth order in the star's 
angular velocity, and the corresponding modes of rotating isentropic stars are gener- 
ically hybrids whose spherical limits are mixtures of axial and polar perturbations. 

These hybrid rotational modes have already been found analytically for the 
uniform-density Maclaurin spheroids by Lindblom and Ipser^^^ and numerically for 
slowly rotating n = 1 polytropes by Lockitch and Friedman and Yoshida and 
Lee^"^-* have obtained these Newtonian polytropic modes to the next nonvanishing 
order in f2. Because the restoring force for these modes is rotation, as it is for 
the axial r-modes, Ipser and Lindblom call them rotational modes or generalized 
r-modes. 

The r-modes of rotating relativistic stars were studied for the first time only 
recently (Andersson ; Kojima^^^; Beyer and Kokkotas^^; Kojima and Hosonuma 
^^)). As in the newtonian spherical isentropic relativistic star has a large 

degenerate subspace of zero-frequency modes consisting of the axial-parity r-modes 

and the polar-parity g-modes. Although isentropic newtonian stars retain a 
vestigial set of purely axial modes (those having I = m), rotating relativistic stars 
of this type have no pure r-modes no modes whose limit for a spherical star is 



*' An appendix in this paper incorrectly claims that no I = m r-modes exist, based on an 
incorrect assumption about their radial behavior. 
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purely axial. Instead, the newtonian r-modes with I = m > 2 acquire relativistic 
corrections with both axial and polar parity to become discrete hybrid modes of the 
corresponding relativistic models. 

In the slow-motion approximation in which they have so far been studied, r — 
modes of nonisentropic stars have, remarkably, a continuous spectrum. Kojima 
shows that the axial modes are described by a single, second-order ODE for the 
modes' radial behavior; he argues that the continuous spectrum is implied by the 
vanishing of the coefficient of the highest derivative term of this equation at some 
value of the radial coordinate, and Beyer and Kokkotas make the claim precise. As 
the latter authors point out, the continuous spectrum they find may be an artifact 
of the vanishing of the imaginary part of the frequency in the slow rotation limit. 
(Or, more broadly, it may be an artifact of the slow rotation approximation.) 

In addition, Kojima and Hosonuma have studied the mixing of axial and po- 
lar perturbations to order f2'^ in rotating relativistic stars, again finding a continuous 
mode spectrum. Their calculation uses the Cowling approximation (which ignores all 
metric perturbations) and assumes an ordering of the perturbation variables in pow- 
ers of f2 which forbids the mixing of axial and polar terms at zeroth order. Again the 
continuous spectrum may be an artifact of an approximation that enforces a purely 
real frequency. 

Finally, several authors have looked at the gravitational-wave driven r-mode 
instability in white dwarfs. Andersson ct al'^^ suggested that the instability may 
limit the spin of accreting dwarfs; and Hiscock^^^ found that an r-mode instability 
of dwarfs conforming to the Wagoner scenario - radiating in gravitational waves the 
angular momentum it gains in accretion - was potentially detectable by space-based 
detectors, with a preference for a shorter baseline instrument like LISA. Lindblom 

however, finds that the growth time of the instability is too slow to allow the 
amplitude to reach saturation in less than 6 x 10^ yr; and he finds a very restricted 
range of mass and accretion rates for which the star could be hot enough for long 
enough that the r-mode could grow to a dynamically significant amplitude. Although 
Lindblom's analysis allows instability in accreting dwarfs larger than 0.9 Mq that 
are spun up by accretion to about 0.85i7x, observable consequences of the instability 
do not seem likely. 

Summary 

In a newborn, rapidly rotating neutron star, an r-mode instability appears to grow 
until it is in the nonlinear regime. It then steadily radiates the angular momentum 
of the star until the temperature falls below the superfiuid transition temperature 
of about 10^ K. By this time, the star's angular velocity will have decreased to 
less than 1/5 the maximum (Keplerian) i7. 

The formation of millisecond pulsars from accretion- induced collapse of white dwarfs 
is apparently ruled out for periods less than about 3 ms. 

The observation of Marshall et al apparently implies that ordinary supernovae pro- 
duce neutron stars rotating fast enough to have unstable axial modes. 
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LIGO II may be able to detect the gravitational radiation emitted by this instability 
in supernovae within 4 to 20 Mpc. 

All of these conclusions are based on linear perturbation theory with no mag- 
netic fields, and work is in progress to decide whether the maximum mode amplitude 
is large enough to justify the dramatic predictions. 
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